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Abstract. We investigate the behavior of non-Euclidean plates with constant 

negative Gaussian curvature using the Foppl-von Karman reduced theory of elasticity. 

Motivated by recent experimental results, we focus on annuli with a periodic profile. 

\^ . We prove rigorous upper and lower bounds for the elastic energy that scales like the 

^^ I thickness squared. In particular we show that are only two types of global minimizers - 

(-H ' deformations that remain fiat and saddle shaped deformations with isolated regions of 

stretching near the edge of the annulus. We also show that there exist local minimizers 
with a periodic profile that have additional boundary layers near their lines of infiection. 
These additional boundary layers are a new phenomenon in thin elastic sheets and are 
necessary to regularize jump discontinuities in the azimuthal curvature across lines of 
inflection. We rigorously derive scaling laws for the width of these boundary layers as 
Q>^ , a function of the thickness of the sheet. 
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1. Introduction 

1.1. Non- Euclidean model of plates 

Laterally swelling and shrinking thin elastic sheets are ubiquitous in nature and industry 
and are capable of forming complex surfaces of various geometries. The examples shown 
in figure [1] include shapes formed by differential growth, thermal expansion, and the 
inhomogeneous swelling of hydrogels. Hydrogels in particular have received a lot of 
attention because of their ability to simulate biological growth by differentially swelling 
when activated by external stimuli such as light p^, solvents [2], warm water [3], and pH 
[1] . The morphology of these structures is a result of the sheet buckling to relieve the 
residual stress caused by the material's resistance to cavitation and interpenetration of 
the material [5j. 

One model of such swelling bodies hypothesizes that the equilibrium configuration 
is the minimum of a "non- Euclidean" free energy functional E^d '■ W^^'^(^3d,IK^) — ^ K 
that measures strains from a fixed three-dimensional Riemannian metric gsa defined on a 
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(a) 



(b) 




(c) 



Figure 1. (a) A hybrid species of the Echeveria plant. The ripphng near the 



boundary is caused by local differential growth and has been the focus of several 
works [3 [71 IHl IS] ■ |(b)| The heating and subsequent drying of potato chips generates 
shapes with a hyperbolic geometry. (c) N-isopropylacrylamide (NIPA) hydrogel 
disk that has undergone controlled shrinking [5^. Hydrogels replicate many of the 
characteristic features present in differential growth and are excellent tools for studying 
such processes quantitatively. 



simply connected domain V-^^ C M^ [SI |TOl [111112]. Laterally swelling thin elastic sheets 
of thickness t can be modelled in this framework by a two dimensional Riemannian 
metric g defined on the mid-surface of the sheet P C M^. That is, V^^ can be 
decomposed as V^d = T> x (— t/2,t/2) and in an appropriate coordinate system gsa 
is given by 



iSd 



The non-Euclidean model of elasticity has grown out of an experiment performed by 
Sharon et. al. that studied the self-similar rippling patterns observed along the edges 
of torn elastic sheets and leaves [13]. This rippling pattern was studied in the strip 
geometry by Marder et. al. [6] and Audololy and Boudaoud [3 |8] by using a metric 
that is localized to one edge of the strip. In particular, in [8] it was shown numerically 
that the self similar patterns could be explained as approximate isometric immersions of 
g with the particular isometric immersion selected as the minimum of a bending energy 
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functional. 

1.2. Reduced theories 

In many applications the thickness of such sheets is much smaller then the diameter of 
V. Consequently, there is considerable interest in obtaining reduced energy functionals 
Et defined on a suitable space of mappings of the mid-surface V into M^ such that 
minimizers of Et approximate minimizers E^d in an appropriate asymptotic limit of 
vanishing thickness. For the classical theory of plates, i.e. when g is the Euclidean 
metric, a family of such reduced theories have been derived through the technique of 
F-convergence [H]. For non- Euclidean plates, there are two reduced theories that have 
recently been rigorously obtained as F-limits by Lewicka and Pakzad [I2] and Lewicka, 
Pakzad, and Mahadevan [15] . 



1.2.1. Kirchhoff model: The first of these theories, which is sometimes called the 
Kirchhoff model, states that if the set Aki of finite bending energy isometric immersions 
is not empty, that is Aki = {x e W^''^{V, M^) : (i^'x)^ . Dx = g} ^ 0, then 

F - lim ^Esd = Eku 
t^o t^ 

with the curvature functional Eki ■ PF^'^(P,]R^) — )• M defined by 

Y 



EKi[x]= <^ 24(l + z/)y^ 



'"' -2K 



dAg if X G Aki 



1-u 

oo if X ^ Aki 

where Y and z/ are the Young's modulus and Poisson ratio of the material respectively, 
H and K mean and Gaussian curvatures of the surface x(D) respectively, and dAg the 
area form induced by g [I2]. This reduced theory captures the intuition that in the 
vanishing thickness limit the mid-surface should deform into an isometric immersion 
with a low amount of bending energy. 

1.2.2. Foppl-von Kdrmdn Model: The second reduced theory, which has been called the 
small slope approximation or Foppl-von Kdrmdn (FvK) ansatz, has also been rigorously 
derived as a F-limit when the following assumptions are met: 

(i) The metric g satisfies the following scaling g = go + t^gi, where go is the Euclidean 
metric. 

(ii) The deformation of the mid-surface x : P — t- M^ satisfies an "in-plane" and "out- 
of-plane" decomposition of the form 

X = i + ti± o ri + t^io Xi 

where x ^ W^^'^(^,M^), rj e PF^'^(T',M), i : M^ -)■ M^ is the standard immersion 
and i_L maps into the orthogonal compliment of 2(]R^). 
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In this context it is natural to define the FvK admissible set by ^ = W^''^{T^,M?) x 
iy^'^(D,]R). When we write that a deformation x : I? — t- M'^ satisfies x G ^ we in 
actuality mean that there exists {x, rj) & A such that x = i + ti± o t] + t'^i o x- The 
in-plane strain tensor 'y : A -^ M^^^ is defined in this approximation by 

7(x, V) = {Dxf + Dx+ {Dvf ■ Dt] - gi, (1) 

which measures to O(t^) the deviation of x from being an isometric immersion. 

The reduced energy under the above assumptions is given by the following F-limit: 

with -EfvK : ^ — )■ M defined by 

where Q : M^^^ — )■ M is the quadratic form 

Q{A) = (1 - u)-^ tr(yl)2 - 2 det(A), 

D^rj the Hessian matrix of second derivatives of t] and dA the Euclidean area form [T5] . 
This energy can be interpreted as approximating E^d by an energy that is the additive 
decomposition of a stretching energy functional S'fvK : ^ — )■ M measuring the lowest 
order elastic energy of the in-plane components of the in-plane strain (Dx)^-Z)x— g and a 
bending energy functional -BfvK : v4 — )■ M measuring the lowest order approximation 
of Eyii. Additionally, in this approximation, the terms Arj and det{D'^ri) correspond 
to the mean and Gaussian curvatures respectfully. The Euler-Lagrange equations 
generated by the variation of -EfvK are a modified version of the classical FvK equations 
and have been used to model morphogenesis in soft tissue P [16] . 

1.3. Motivation for and physical implications of this study 

There have been several theoretical ^7\ [IHl |19] and experimental [20] investigations 
that have studied non-Euclidean plates when D is a disk or an annulus with a two 
dimensional metric gxo that has corresponding constant negative Gaussian curvature 
Kq. Since ^xi 7^ for these metrics and domains [211 [22] . that is there exist exact 
isometric immersions with finite bending energy, it follows from the Kirchhoff model that 
in the vanishing thickness limit the mid-surface of the minimizers of E^^ should converge 
to a fixed shape. Furthermore, in a related work Trejo, Ben Amar, and Miiller [23j have 
shown that surfaces with constant negative Gaussian curvature can be constructed from 
closed curves generating pseudospherical surfaces which has implications to the modeling 
of growing bodies with a tubular topology. 

Experimentally, upon swelling, annular hydrogels with a fixed inner and outer radius 
and a prescribed two dimensional metric gxo obtain a periodic profile with n waves. The 
number of waves refines with decreasing thickness according to the scaling n ~ t~^/^ with 
the bending energy ii^xi diverging according to the scaling E^i ~ t~^ [20]. Therefore, at 
least for the range of thicknesses considered in the experiment it appears that the annuli 
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are not converging to a fixed shape. If we believe tliat tlie experimental configurations 
are global minimizers of the elastic energy, then they should converge weakly to an 
isometric immersion with finite bending energy. It is this apparent contradiction between 
the experimental observations and the theoretical prediction that motivates the study 
in this paper. 

A possible explanation for this apparent disagreement between theory and 
experiment is that the range of thicknesses which are accessible to the experiments 
corresponds to an intermediate asymptotic regime, different from the regime of vanishing 
thickness. That is, although the number of waves increases in the intermediate regime, 
for smaller thicknesses the system cross over into a different regime at which point the 
shape (and hence also the number of waves) will stabilize and converge to an exact 
isometric immersion. 

To study this question, we need to determine the lowest energy configurations of 
the sheet not only in the limit t — )• 0, but also for alH > 0. In addition to the thickness, 
a second relevant parameter in this system is the dimensionless curvature of the sheet 
e, namely the curvature of the sheet normalized by its radius. For the experiments, 
this parameter is smaller than 1. Thus it is appropriate to study these sheets in the 
Foppl-von Karman (geometrically linear) regime |15] . 

We obtain rigorous scaling laws for the Foppl-von Karman elastic energy for all t > 
and all configurations with n- waves. In particular, since we are interested in all t > 0, 
we do not assume that the configurations are close to isometric immersions. Our scaling 
laws for the energy are ansatz-free, and applicable to all n-wave configurations. The 
rigorous bounds show that the n = 2 configurations have lower energy than all the n > 2 
configurations for all thicknesses t > 0. This falsifies the idea that the experimentally 
observed n > 2 configurations are global minimizers for the elastic energy. 

In addition to the global minimizers (n = 2) we also determine the local minimizers 
{n > 2) of the elastic energy for t > 0. As t — )• 0, these local minimizers also converge 
to isometric immersions which however are not smooth. These isometric immersions 
do not self-intersect and for this metric there is no analogue of the skewed e-cones 
introduced and studied in [23]. But, all the local minimizers with n > 2 waves converge 
to configurations with singularities ( "lines of inflection" ) corresponding to a continuous 
tangent plane, but a jump in the curvature [I9]. This should be contrasted with the 
elastic ridge singularity in crumpled sheets which mediates a jump in the tangent plane 
[25]. 

We have determined the shape of the local minimizers for t > by solving for the 
boundary layers which resolve these jumps in curvature. For n > 2, we identify two 
types of overlapping boundary layers corresponding to reductions in the Gauss curvature 
and the mean curvature respectively. This overlapping is, to our knowledge a unique 
phenomenon and one that is potentially observable in experiments. 

Finally, we observe that although the n > 2 local minimizers are periodic, the 
fact that they have singularities (or narrow boundary layers) implies that they do not 
have a sparse representation in terms of a Fourier basis. Thus, these structures will be 
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difficult to analyze numerically using a Fourier basis expansion as was done for the strip 
geometry by Audoly and Boudaoud [8]. 

1.4- Structure of paper and main results 

In this section we describe our results in a mathematically precise manner and also 
outline the structure of the paper. As discussed above, we study the convergence and 
scaling of local and global minimizers of the Foppl-von Karman energy with the metric 
gKo for decreasing values of t > 0. A major point of this paper is that we prove that in 
the FvK approximation that for all thicknesses saddle shapes {n = 2) are energetically 
preferred over deformations with more waves (n > 2). 

Motivated by the experimental results in [20] we focus on studying minimizers 
with a periodic profile when V is an annulus with inner and outer radius po and R 
respectively. To be precise, for n G {2, 3 . . .}, we study minimizers over the admissible 
set of n-periodic deformations An C A defined by (x, v) ^ -^n if and only if in polar 
coordinates (p, 9) the out-of-plane displacement rj satisfies 

(i) rj is periodic in 6 with period 27r/ra, 
(ii) 7] vanishes along the lines 6 = and 6 = ii/n, 
(iii) r]{e-TT/n) = -r]{e). 

We call the lines 6 = mn/n, m & {0, . . . ,2n — 1}, lines of inflection. 

In section 2 we construct the FvK elastic energy by expanding g^o in the 
dimensionless curvature e defined by 

e = ^/^K'oR. (3) 

We repose the variational as the minimum of a dimensionless functional Sr = S + t^B, 
with the dimensionless thickness r defined by 

r = t/(V3i?e). (4) 

We further derive the Euler-Lagrange equations and natural boundary conditions 
corresponding to the variation of Sr- 

In sections 3 and 4 we construct minimizers over the admissible sets Af = {x G 
A : i3[x] = 0} and Aq = {^ E A : 5[x] = 0}. The set Af consists of flat deformations 
and the minimum over Af is obtained by solving the Euler-Lagrange equations and its 
natural boundary conditions with rj = 0. The set ^o consists of isometric immersions 
in the Foppl-von Karman approximation and can be obtained by solving the Monge- 
Ampere equation 

det{D^7]) = -1. (5) 

The global minimum of B over solutions to (j5]) are obtained by the approximate minimal 
surfaces rj = xy and thus following the results in [T9] we have the following theorem: 

Theorem 1.1. Let x.* E A be a sequence with corresponding out-of-plane displacement 
rj* such that inf xg^fi £^t [x] = ^r[x*]. Then, 



Defects and boundary layers in non-Euclidean plates 7 

(i) lim inf 4^r[x] = lim 4^r[x!] = inf B\x\ = 27r(l - plR-^). 

('zzj There exists a subsequence rj*^ and x* G >4.o wt/i corresponding out-of-plane 
displacement rj* such that rj*^ -^ rj* . Moreover, there exists A G 5*0(2) and 6 G M 
such that ri*{A{x, y)) + b = xy. 

This theorem proves that in the vanishing thickness hmit saddle shaped deformations 
are energetically preferred. 

Following the construction in [T9|, we show that An fl ^o 7^ 0, that is there exists 
n-periodic isometric immersions, and using these constructions as test functions we show 
the following result: 

Lemma 1.2. Let x* G An such that £t[^*] = iiifxeyl„ ^r[x], then 

Sr[^*] <mm{Cn^T\T}, 

where T = infxe^^ ^r[x] and C is a constant independent of n and t. 

This lemma captures the difference in scaling of £r with r for fiat deformations and 
isometric immersions by using specific test functions on Sr- Furthermore, this lemma 
illustrates the growth in n of the total energy for n-periodic isometric immersions. 

In section 5 we numerically minimize Sr over An using a Rayleigh-Ritz type method. 
The results of the numerics indicate that minimizers transition from fiat deformations 
to buclded shapes such that with decreasing thiclcness the deformations converge to 
X G Anf^Ao. Moreover, the results of the numerics indicate that for all thickness values 
the 2-wave saddle shape is energetically preferred over profiles with a higher number 
of waves. In the bulk of the domain the stretching energy of these buckled shapes 
is approximately zero with regions near the edges and, for n > 3, along the lines of 
inflection in which stretching energy is concentrated. These results indicate that with 
decreasing thickness the stretching energy is concentrated in shrinking boundary layers 
in which the bending energy of the isometric immersion is reduced by adding a small 
localized amount of stretching. 

In section 6 we improve on lemma [L2l and rigorously prove ansatz free lower bounds 
on Sr- Specifically, we prove the following theorem: 

Theorem 1.3. Suppose n G {2,3,...} and r > 0. There exists constants c and C 
independent of n such that 

min < cnT"^, — > < inf S\x\ < minjCn^r^, J-"). 
[ 2 J ~ ^eAn ^ ' ~ ^ 

This theorem confirms the numerical results that with decreasing thickness the elastic 
energy of the minimizers scales like that of an isometric immersion. Furthermore, since 
both the upper and lower bounds in the theorem grow with n it quantifies how the elastic 
energy is penalized by adding more waves to a deformation. Additionally, this theorem 
proves that for all thicknesses the 2-wave saddle shape is energetically preferred. 

In section 7 we determine the scaling with r of the width of the boundary layers in 
which the stretching energy is significantly non-zero. Near the edges of the annulus the 
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out of plane displacement can be additively decomposed into two terms that separately 
lower contributions to the bending energy from the mean and Gaussian curvatures. The 
width of the regions in which these terms are relevant satisfies the following scaling: 

(i) The boundary layer in which one term lowers the magnitude of the Gaussian 
curvature has the following scaling 

width{6) p=pg^R ~ t2|i^Q|-V4pg(, / _\ cos{9) cos{n/n - 6). 

(ii) For n > 3, the boundary layer in which one term lowers the magnitude of the mean 
curvature has the following scaling 



width(^)p=PQ_ij ~ ^H-^ol ^^"^ysin ( — j sec(6') sec ( 



The overlap of these two regions forms the complete boundary layer in which the total 
bending energy is lowered by allowing localized stretching. The width of these boundary 
layers has the same scaling in thickness obtained by Lamb [26j in the context of vibrating 
shells and Efrati et. al. for non- Euclidean plates with Kq > [H]. The fact that there 
is no boundary layer to reduce the mean curvature in the case n = 2 is a consequence 
of the fact that rj = xy is a. minimal surface in the FvK approximation. 
The width of the boundary layers near the lines of inflection scale like 

width(p)^=o ~ ^^P^I-^oT^; 

which has the same scaling with thickness, but not p, for minimal ridges formed 
by crumpling [25]. In this boundary layer the mean curvature is locally reduced by 
correcting a jump discontinuity in the azimuthal curvature. This reduction of energy 
near this type of singularity is different from the regularization near a ridge singularity in 
which the bending energy diverges while the stretching converges to zero with decreasing 
thickness [271128]. 

We conclude the paper in section 8 with a discussion of the implications of the 
results of this paper and future mathematical directions stemming from this work. 

2. FoppI - von Karman equations for metrics with constant negative 
Gaussian curvature 

2.1. Fdppl - von Karman ansatz and elastic energy 

Let ^"xg denote the hyperbolic plane with Gaussian curvature Kq and let V C Mkq 
be an annulus of inner radius po and outer radius R centered at the origin of M^o- 
Introducing the dimensionless curvature e defined by 



\Ko\R (6) 

it follows that in geodesic polar coordinates (p, 6), the metric on V is [29] : 

, 2 sinh^ (ep/i?) ,„2 
g = dp2 + , ^// ' de^. 
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Therefore, expanding in e, g can be approximated to order e^ by g = go + e^gi where 
Qo is the Euchdean metric on M^ and 

^^ ^ 3R^ ^^^ ^ W^ ^^^ ^^^ ~ ^^^ ^"^^ ^ "^ ^^^'^ ' 
with the "Cartesian" coordinates u = pcos{6) and v = psin{6). 

A deformation x : 2) — )■ M'^ can be constructed that is an isometric immersion of g 
to order 0(e) by assuming that 

X = i + ei± o 7] + e'^i o X, (7) 

where x ^ W^'^{V,M.'^) and rj G iy^'^(T',M). In this context, when we write that a 
deformation x : P — ;■ M^ satisfies x G ^ we mean that there exists (x, f]) & A such that 
x = i + ei± o 7] -\- e^i o x- Based on ([2]) we define the elastic energy i^^ : ^ — ;■ M by 

Ye' ' 



E,h] 



8(1 + i^) Jv 



tr(7)^-2det(7) 



1-z/ 

dudv, . 



t^ / 1 

■tT{D^T]y -2det{D^ri) 



where again 7 is the in-plane strain defined by ([T]) and D^rj denotes the Hessian matrix 
of second partial derivatives of 77. A deformation x G ^ is called an equilibrium 
configuration in the FvK ansatz if 

Et[^] = MEtly]. 

Define the dimensionless variables x,y,x', and rj' by 

u = Rx, V = Ry, r = Rp, x = R'X' 1 V = Rv' ^ 

and (-D^)' the Hessian operator in the coordinates x and y. Then {D^Yr]' = D^rj and 
thus introducing the dimensionless thickness r defined by 

T = t/ (y/?>R^ (8) 

and dropping the / notation we have that 



tr(7)2-2det(7) ) dxdy 



Ye'R^ 

1 



+ r^ f f_L_tr(^{D^yv')-2det ((D^)'^)) dxdy. (9) 

where B is an annulus with inner radius ro = po/R and unit outer radius. Minimizers 
of Et also minimize Sr and therefore we will restrict our attention to the functional Sr- 
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2.2. Euler- Lagrange equations 

Throughout the rest of the paper let dB denote the boundary of B with outward normal 
vector field n and tangent vector field t. In this subsection we will assume that x G ^ 
extremizes £r with out-of-plane displacement rj and in-plane displacement x = (Xi) X2)- 
Taking the variation of Sr with respect to Xi "we obtain: 



IB 



ih) ^^'' + ^''^ ^^^^ ' ^ (^^^^^^^ ' ^'Wy^''' 



dxd|/ 



= 4 / / ^1_ — ^^ ^^2 j ■ n^xi ds - 4 / V • ( ^y— — ^,712 j Sxi dxdy. 
Similarly, the variation with respect to X2 yields: 

5x2-^r[x] =4 / (712, ^^^_^" ' j ■n(5x2ds-4 / V- f 7i2, ^yTT^' ) ^Xi dxdy. 
Therefore, the governing equations for the in-plane strain are 

V-(2^i^Y^,7i2)=0, V ■ (712, ^^2^^^) = 0, (10) 

with the natural boundary conditions 

„ /^ 711+1^722 ^, \\ _n T. f'o; i^7ll+722 N I _n 

n-(^ i„^ '^i2J|g^-0, n-(^7i2, ^_^ JUij-U- 

Furthermore, (TTOl) implies that there exists a potential function $ G W^^'^(i?,]R) 
satisfying 

jj2^^J_f ^711+722 -(l-^)7l2y .^^. 

1 - Z/ y -(1 - Z/)7i2 7ll + 1^722 J 

Now, if we assume 7 is second differentiable, then differentiating and eliminating 
terms involving x from ([1]) we have the following compatibility condition between 7 and 
r/: 

9 7ii <9 7l2 f^ 722 

(9y2 (9x9?/ (9x2 
where the operator [-,■]: C2(fi,R) x C^{B,R) ^ M is defined by 

2 \ dx"^ dy"^ dy"^ dx"^ dxdy dxdy J 

Therefore, inverting (ITT]) and substituting into (TT2|) we have the following proposition. 

Proposition 2.1. Ifx* G >4. wi/i potential $ anc? out-of-plane displacement rj , Sr[x*] = 
inf xe^ £^r [x] , anc? z/ $ zs (weakly) four times differentiable then $ satisfies the first 
Foppl - von Kdrradn equation: 

^ A2$ + [r/,r/] = -l. (13) 



2(1 + z/ 
Also, <l> satisfies the following natural boundary conditions: 
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Furthermore, by the definition of the potential $ given in by f llip the stretching 
energy can be expressed as: 

1 



5ixl 



(A$)^_2[$,$] dxdy. 



'B x^^i' 

The next proposition allows us to simplify this energy further and will be useful in the 
next section. 



Proposition 2.2. If x* E A with potential $ and infxe^i^rM = ^r[x*] then 



5[x1 



(A$)^ dxdy. 



Proof. Let x* E A with corresponding potential function $ such that infxg^^r[x] = 
£^T-[x*]. Using Stokes' theorem and applying the boundary conditions (IT^ . we have that 

dx"^ dy"^ dxdy dxdy J 



[$, $] dxdy 



B 

9x \ 9x dy"^ dy dxdy J dy \dy dx^ dy dxdy J 

— — 1 • n drS 

QB \ dx dy"^ dy dxdy ' dy dx"^ dx dxdy J 



dxdy 



QB dx \ dy"^ ' dxdy 






■ nds 



Therefore, 



5[x*] = / {-^ {^^f - 2[$, $]) dxdy = / -^ (A$)2 dxdy. 



D 



Now, taking the variation of S with respect to rj and applying the boundary 
conditions (ITU) we have 



(5^5 [x] 



B 



1 , , f dri d dr] d 

1 — u \dx dx dy dy 

dri d ^ dr) d ^ dri d ^ dr) d ^ 

dy dy dx dx dx dy dy dx 



dxdy 



9^$ dr] 9^$ dr] 9^$ dr] 9^$ dr] 



B V dy'^ dx dxdy dy ' 9x^ dy dxdy dx 

dr] /(9^$ 9^$ \ dr] f' 9^$ 9^$ \ 

9x \ dy"^ ' (9x9y / dy \ dxdy ' (9x9y / 



dB 

4 



V(5r7 dxdy 

^r^ds 



[$,r7]5?7dxdy 



B 



4 / [<l>,r7]5?7dxdy. 
'b 
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the variation of B with respect to t] is 



Furthermore, following the derivation in 
given by 

1 



S,B[^] 



A rjSri dxdy + 

B JdB 

1 dAr] d 



1 



At] — n ■ D T] ■ n 



d6r] 
dn 



ds 



as 



+ (n^.Z)2^.t) 
1 — z/ an at ^ 



5rjds, 



where is the angle n makes with the x-axis. Therefore, we have the following 
proposition. 

Proposition 2.3. //x* G A with potential ^ and out-of-plane displacement rj , Sr[x*] = 
infxe^^r[x], and if rj is (weakly) four times differentiable then 



[$,r/] 



-A^r/. 



(15) 



4(1 -z/) 

This equation is called the second Foppl - von Kdrmdn equation. Furthermore, rj 
satisfies the following natural boundary conditions: 

1 



l-v 
1 dAri 



Arj — n ■ D T] ■ n 



0, 



dB 



d 



+ Tsr (n^ • D^V ■ t) 
1 — z/ an at 



0. 



(16) 



:i7) 



dB 



Remark 2.4. The equations [W\) and [TB) along with the boundary conditions ( [7^ , [TB) 
and ( [77| j are a complete system of equations governing the deformation of the sheet. For 
general geometries, the boundary conditions are intractable and in later sections we will 
use the geometry of B to simplify these equations further. 



3. Flat solution 

Define the admissible set of fiat configurations Af C .4 by Af = {x G ^ : r^ = 0} 
and define F = infxg^^ ^r[x]. If x G Af extremizes Sr with potential function $ then, 
assuming radial symmetry, we have by the first FvK equation (fT3|) and the boundary 
condition JHM that $ satisfies the following boundary value problem 



+ - 



ar^ 



/'i/y "^ /Y"^ /iT' 



1 52$ 1 5$ 

+ ^^^ = -2 1 + z/ 



'lO /if^ 



dr 



0. 



(18) 



r=rnA 



The general solution to this differential equation is 

4 
$ = Ci + C2r^ + C3 ln(r) + c^r"^ ln(r) - (1 + '^)^, 

where Ci, 02,03 and C4 are arbitrary constants. Since Ci and C3 ln(r) are harmonic 
functions it follows from proposition 12.21 that we can assume without loss of generality 
that Ci = C3 = 0. Consequently, the solution to (!T8|) is 
l + iy 



$ 



321n(ro) 
Therefore, we have the following result: 



[(r^ - 1 - 2 In(ro)) r^ + 2 (r^ - l) r^ ln(r) - In (ro) r^] . 
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Lemma 3.1. Suppose x* G Af satisfies £t[^*] = infxey^f ^rM- Then, 



£r[^*] = J^ 



o W - 1) + {rl - 1) W + 1) In(ro) + - (l - r',) In(ro)^ 



321n(ro)2 [2 

4. Isometric immersions 

Define the set of isometric immersions by ^o = {x G ^ : 7 = 0}. By tlie 
compatibility condition ( IT2l) . it follows that the solvability condition for the equation 
7 = is the Monge- Ampere equation 

[r],r]] = det{D^=-l, (19) 

which is a version of Gauss's Theorema Egregium in the FvK ansatz. Consequently, if 
X G ^0 with corresponding out-of-plane displacement rj and x satisfies flT3l) and flTSl) 
then A'^rj = 0. Therefore, it is natural to look for solutions to the FvK equations by 
finding solutions to the above Monge-Ampere equation (fT9!) satisfying A'^rj = 0. 

4-1- Saddle isometric immersions and convergence of minimizers 



By the definition of Sr (Ej) and the Monge-Ampere equation (1191) it follows that for all 
X G ^0 

^ inf SJ-k] = inf [ -^ {Arif dxdy + 27r(l - r^). (20) 

T^ xe^o xe^o J^ 1 — z/ 

Consequently, since tj = xy satisfies (IT9l) and is harmonic, it follows that the global 
minimum of f l20|) over ^0 is obtained. Therefore by the arguments presented in [19] we 
have the following results concerning the convergence and scaling of minimizers of £r'- 

Theorem 4.1. Let ^*^ E A he a sequence with corresponding out-of-plane displacement 
rj* such that infxg^<£^rM = ^r[x*]. Then, 

T->-o xe.4 r^ T-s>o r^ xeAi 

(a) There exists a subsequence rj*^ and x* G ^0 with corresponding out-of-plane 
displacement rj* such that t]*^ -^ rj* . Moreover, there exists A G 5*0(2) and h eM. 
such that rj*{A{x, y)) + h = xy. 

Proposition 4.2. Let x* G ^ such that £t[^*] = infxe^^r[x]. Then, 
£r[x*] < min{27rr2(l - r^), jr}. 

Remark 4.3. There are two natural test functions for the elastic energy, flat 
deformations and isometric immersions, the upper bounds in the above lemma 
correspond to the minimum of the elastic energy over these two classes of deformations. 
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4-2. Periodic isometric immersions 



Let 



^n = 1/ U^ - cot 



y 



CSC 



COS 



29 



COS 



(21) 



which satisfiesthe Monge- Ampere equation (fT9l) . A^r/ = and ^„ = along the hnes 
6* = and 6 = n/n. Consequently, by taking odd periodic extensions of ?]„ (see figure [2]) 
an ra-periodic isometric immersion rjn can be constructed [19]. Therefore, An H Aq y^ 
and by calculating the energy of this configuration we have the following result which 
agrees with proposition 14.21 when n = 2. 









Figure 2. 1. The one parameter family of isometric immersions with out-of-plane 
displacement rj = y {x — cot (tt/u)) parametrized by n e {2, 3, . . .} vanishes along the 
lines 6 = and 9 = Tr/n. 2. We can "cut out" the section of the surface bounded 
between these two lines. 3. We can take the odd reflection of this isolated piece of 
the surface about the the line 9 = 0. 4-7. By continually taking the odd reflection of 
sectors about lines where the surface vanishes we can construct a periodic isometric 
immersion. 



Lemma 4.4. Let x* G An such that Srl^*] = iiifxey4„ <S't[x], then 

2 ^'47rcot^(7r/r2) 



•^rfx*] < min < r 



where C is a constant independent of n and t. 

Proof. The middle term in the chain of inequalities follows from calculating the bending 
energy of the deformation x„. Expanding near n = oo we have that 



47rcot2 (^) 



+ 27r 



4n^ 



+ 27r 



in 



+ ... . 



1 — v ' 7r(l — z/) ' \y 3(1 — ntt) 

Consequently there exists M > and Ki > such that if n > M then 
47rcot2(^) 



1-u 



+ 2tt < Kin^. 
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Differentiating we have that 

dn \nJ n^ 

and thus cot^ ( — ) is monotone increasing. Consequently, on the interval [2, M] we have 
that 

471 

~l-u ' — - 1 i_j^ '2 
Therefore if we set C = ina.x{Ki, K2} then the result follows. D 

5. Numerical solutions 

The boundary value problem given by the FvK equations and its natural boundary 
conditions (11311171) can be approximately solved by numerically minimizing Sr- To do 
this, it is convenient to write the energy in the form 



47rcot2(^) /vrcot^f^) tt , , . 



' B 



+t' 



{ill + 722) + 711 + 27L + 722 



dxdy 



V 



(Arif + \D^ri\ 



dxdy. (22) 



'B 

Sr can then be discretized using a finite difference scheme to approximate the integrand 
and a quadrature rule to approximate the integrals. This discretization of fl2^ generates 
a sum of quadratic terms that can then be minimized using Matlab's minimization 
routine Isqnonlin PT] . 

In figure [3] we plot the elastic energy of numerical minimizers of Sr for decreasing 
values of r with vq = 1 and z/ = 1/2. The data in this figure was generated using the 
algorithm outlined in the previous paragraph with a 40 x 40 mesh in the coordinates r 
and 6, a fourth order approximation of derivatives, and a second order approximation 
of integrals. The solid line corresponds to the upper bound in lemma 13.11 for flat 
deformations while the thin dotted line is the upper bound in proposition 14.21 for 
isometric immersions. In this figure three representative minimizers are plotted and 
coloured by the Gaussian curvature K = det{D'^r)) for various values of r. These 
three surfaces and the scaling of the energy illustrate that with decreasing thickness 
the minimizing surface transitions from being fiat to one that is close to the isometric 
immersion rj = xy. For the left most surface, the regions in which the Gaussian curvature 
is substantially different from —1 are localized to the edges of the annulus and shrink 
with decreasing thickness. This indicates that with decreasing thickness the stretching 
energy is being concentrated in boundary layers in which the bending energy of the 
isometric energy is slightly reduced. 

In figure H] we again plot the numerical minimizers of Sr using the same parameters 
and discretization used to generate figure [3] but with the boundary conditions 77 = 
along the lines 6 = 0, vr/ri, 27r/n, . . . for n G {2,3,4,5}. These boundary conditions 
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Figure 3. The normalized elastic energy of numerical minimizers of £r with ly — 1/2 
and To = 10""'^ for decreasing values of r. The solid line corresponds to the elastic 
energy of minimizers over the set of flat deformations given by lemma [33] The dotted 
horizonal line corresponds to the global minimum of the elastic energy over the set 
of deformations satisfying 7 = 0. The three configurations plotted are colored by the 
approximate Gaussian curvature K = [77,77]. These three surfaces illustrate that with 
decreasing thickness the surface transitions from being a flat surface to one that is 
close to the isometric immersion j] — xy with localized regions of stretching near the 
inner and outer radius of the disk. 



were selected to generate numerical minimizers over An- The dotted horizontal lines 
correspond to the upper bounds for n-periodic isometric immersions given by lemma 1131 
while the solid hne is again the upper bound in lemma 13.11 The four plotted surfaces 
are coloured by the approximate Gaussian curvature K = [77, 77] and were selected to 
compare the geometry of the boundary layers for n > 3 with the case n = 2. For n > 3 
additional boundary layers form around the lines of inflection in which the surface 
stretches to reduce the local mean curvature of the surface. 

6. Scaling laws for the elastic energy of periodic configurations 

In this section we derive ansatz free lower bounds for n-periodic configurations. The 
essential idea of this section is that the bending energy near the edge of the annulus 
controls the stretching energy in the bulk of the domain. Before we state and prove 
this lower bound we need several intermediate results. Let n G {2, 3, . . .} and define the 
following sets: 

(i) Sn is the sector in R^ bounded between 6 = and 6 = ir/n, 

(ii) For z > Tq, B^ G B is the annular region defined hy B^ = {{r,6) E B : Tq < z < 
r<l}, 
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Figure 4. The normalized elastic energy of numerical minimizers of £r over An with 
To = 10~^ and v = 1/2. The dotted horizontal lines correspond to the upper bounds 
given by lemma 14.41 while the solid line corresponds to the upper bound in lemma 
13.11 The configurations plotted are colored by the approximate Gaussian curvature 
[r],ri]. In the vanishing thickness limit minimizers of £r converge to an n-periodic 
isometric immersion. The numerical minimizers have localized regions of stretching 
near the inner and outer radius and along the lines of inflection. The existence of these 
regions indicate that the minimizers are perturbations of an isometric immersion with 
boundary layers to account for the natural boundary conditions. 



(iii) B^ = Sn^i B^ is a wedge shaped region bounded between 9 = 0, 9 = n/n, r = 1, 
and r = z. 

Furthermore, let r* = niax{.95,ro}. 

Lemma 6.1. //x G An with corresponding out-of-plane displacement t] then 



drj 
dr 



dA< 



^ Jb^ 



V drd9 ) 



dA 



(n) 
(iii) 



ri^dA< 



TT 



48^4 ./^, V 89'^ 



d'^Tj^ 



m) -12^2/1^^2 



dA 



dA 



(23) 
(24) 
(25) 



Proof. Since x G An we have that ?7 = along the lines ^ = 0, Ti/n and thus 
drj 



dr 



0. 
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Consequently it follows from Poincare's inequality with the optimal constant [22] that 

ydr J ~ ri^ Jq \drd6 ) 
Therefore, integrating we have that 



Lm^^<-hLm-^ 



proving item (i). 

Constructing the Green's function G{6^ (p) for the operator -^ with Dirichlet 
boundary conditions on [0, vr/n] we have that 

d'^ri 



v{r,e)= / G{e, 



where 



Therefore, 



""^''^^-^ ^A<P-rr/n), if < ^ 



V\r,e) < I / G\eA)dA 1 I / ( ;;;^ 1 d^ 1 = ^\ ^^^' /" ( ^] dO 




3nn In V 99 



rh^A<J^( fSVd^. 



48n3 7o 
Integrating we have that 

proving item (ii). 

Differentiating, we have that 



89 Jo 89 8(t)'- 



rd0. 



Therefore 




f ^^ \ m\ i [" f '^^V \^ -I a\ f '^ n ^^^ ^ /" " /^ ^^^ 



Km}'' \k la^^j^^hU-"^- 



TT 
< 



I2n jQ 




Integrating we have that 



??VdA<X/ (?2] dA 



iB^ \89j - 12n2 Jb. {89"^' 

proving item (iii). D 
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Lemma 6.2. Let x G An with corresponding out-of-plane displacement t] . If B[x\ < Bq 
and n > 2 then there exists a constant C independent of n and x such that 



_B'* 






\drde) 



dA<C 



B. 



n 



Proof. Since i3[x] < Bq it follows that 
D^r^\MA -- 



Therefore, 



1 fldr] d'^ri \ ' 
.2 \j-Q0 drdO J 



d^rj drj 
.4 \gQ2+^-g;: 



Br* 



1 dri d'^T] \ ■ 
r 86 drdO J 



d'^rj drj 

QQ2 Qj. 



dA< 



+ 



So 
2n 



drj 
dr 



dA< 



2n' 



and thus applying the elementary inequality (a + b)^ > \a^ — 26^ we have that 



1 

2 .1 B^" 



d'^rj 



\drde 



dA-2 



Br* 



L(^\ +r^(^ 



09 



dr 



- 2n 



Consequently, applying lemma EH] we have that 



Br* 



d'^T]\ / d"^!] Y 



06^ 



\ drdO I 



dA<^ + 4 

n ./fir* 



,*N2 



<^+ 



TT 



d'^ri 



dA+ — 



n 3n'^(r*Y jRr* \d6^ J n- ;nr 

Therefore, since r* > 7r/(2-\/3) and n > 2 the result follows. 






dA 



D 



Lemma 6.3. Lei x G ^„ wi/i corresponding out-of-plane displacement rj. //i3[x] < Bq 
and n > 2 then there exists a constant C independent of x and n such that 

|Vr/|MA < C-^° 



Proof. By lemmas 16.11 and 16.21 it follows that 

^4 r '' d'^rj^ 



r]^dA < 



IT 



Br 



A8n^ 



Br 



de^ 



dA < c,% 



|Vr/|MA = 

<C, 



d'^r] 
dr 



+ 4f?? 






oe 



dA< 



Br* 



B, 



w 



By rotational symmetry of An we have that 



Br* 



r]'^dA<Ci^ and 



W" 



d'^rj 
dr 



\VT]\^dA < C2 



B, 



Br* 



n^ 



+ 






drj 



dA 



Defects and boundary layers in non-Euclidean plates 20 

and therefore by a multiplicative inequality [33] it follows that there exists a constant 
C3 independent of x and n such that 

/ \Vv\^dA<cJ [ v^dAY ( [ {ri^ + \Vr]\^ + \D^vndAY 

J B^* \Jb^* J \Jb^* / 



B? 



n"^ \ n* n^ 



3 



BS / Bo , ^ So ^ ^ 










D 



Lemma 6.4. Let x G A2 with corresponding out-of-plane displacement 7] . If Br[^] < Bq 
then there exists a constant C independent of x such that 



[ \Vr]\^dA < C 



B"^ 



Proof. Since x G ^2 it follows that the means of the functions r/, g^, and -t^ satisfy 

' ^ vdA = 0, ^^^^[ '^dA = 0, ^-^r^f '^dA = 0. 



7r(l - (r*)2) y^.. ' ' 7r(l - (r*)2) 7^.. dx ' 7r(l - (r*)2) ^,.. 9y 

Therefore, by Poincare's inequality it follows that there exists constants Ci and C2 
independent of x such that 

p-] dA<Ci [ \Vr]\MA and / \Vr]\MA < C2 [ \D^v\^dA. 

Bt* \0X J J Br* J Br* J Br* 

Consequently, by applying the same multiplicative inequality as in the proof of lemma 
16.31 the result follows. D 

Lemma 6.5. Let x G An- If B[x] < Bq then there exists a constant C > independent 
of X and n such that 

s[^]>]-:f-c^. 

Proof. Let 77 be the out-of-plane displacement corresponding to x and x = (/? 9) the 
in-plane displacement with component functions f,g & W^^'2(B, M). Then, by lemmas 
16.31 and 16.41 it follows that there exists a constant C independent of x and n such that 

i2" 



[ iVr^lMA < [ |Vr/|Mv4 < C- 
Jb-* Jb n- 

Therefore, the following inequalities hold: 
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(i) 

7?,dxd,= / I 2^ + f |?V - l^^l d.dy 



B' J B 



dx \dx/ 3 



111 



2 ^ /o\2/o\2 



2 ./fir* V 5x 3 y "^ n^ 



to 

2 J^r* \ oy 6 J ri'^ 



(ii) / 72^2 dxdy > i / f2^ - ^ V dxdy - 2C^° 



^r 



B'- 






B"- 



~ 2 yB-* V ^y dx 3 



\ 2 ^^2 

da;d2/-2C^. 



(iv) 



I [ f df dg x^ y^Y , , ^i32 



hn + 72.)^ dxdy > - J^^^ ^2^ + 2g - ^ - ^J d.d, - 8C^. 
Therefore, since 5[x] can be rewritten in the following form 

5[x] = / (- tr(7)2-2det(7) J dxdy 

= J (yZ^ ^^11 + ^22)' + 7n + 27?2 + 722) dxdi/ 

> f ^ (yZ7 (^11 + ^22)' + 7n + 27?2 + 722) dxdy 
it follows by items i-iv that there exists a constant C independent of n and x such that 



5|x|>l^-4 



n 



Remark 6.6. The preceding lemma is the essential estimate that quantifies the trade off 
between bending and stretching energy. Furthermore, it shows that in the limit n — )■ 00 
that S is bounded away from zero. That is, as the bending energy increases by adding 
more waves there is no reduction in the stretching energy. This is in contrast to the 
behaviour of a minimal ridge in which with decreasing thickness the bending energy 
diverges while the stretching energy converges to zero [2j 
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Theorem 6.7. Suppose n G {2,3,...} and r > 0. There exists constants c, C > 
independent of n such that 

niin{ cnr^, J-"/ 2} < inf £^[x] < niin{Cr2^r^, J-"}. 

X6,4n 

Proof. By lemma 16.51 there exists a constant ci independent of n and x such that 

Therefore, minimizing the function £ = S + t'^B subject to the constraints 

5>0, r2i3>0, 5>i^-cif, 
it follows that there exists a constant c independent of n and x such that 
max I cnr^ , J-'/2 1 < inf £ [xl . 

The upper bound follows from lemma 14.41 D 

Corollary 6.8. There exists n* > 2 such that if t < [J^/(47r(l — Tq))] ' and n > n* 
then 

inf £r\^ < inf £'r[x]. 

Proof. By theorem 16.71 it follows that there exists a constant c independent of n and 
r such that minjcnr^, J^/2} < infxe^„£^rM. Furthermore, if r < [T/{Att{1 — r^))] 
then 27r(l — ro)r^ < J^/2 and if n > 27r(l — rQ)/c then 27r(l — tq)^ < en. Therefore, 
letting n* = 2tt{1 — rl)/c it follows that if r < [J'/(47r(l — Tq))]^'^ and n > n* then by 
proposition 14.21 

inf ^^[x] < 27r(l - tDt"^ < min{J^/2, cnT"^} < inf ^^[x]. 
xeyl xeAi 

By theorem 16.71 and the fact that the upper bound in proposition 14.21 corresponds to 

the elastic energies of deformations in A2 it follows that min{2cr^, J^/2} < 27r(l— ro)^T^. 

Therefore, c < 7r(l — Tq) and consequently n* = 27r(l — rD/c > 2. D 

Corollary 6.9. Let n G {2,3,...}. There exists constants c,C > independent of n 
and T such that if r < J-' / {2c)n^^^'^ then 

criT^ < inf £r[^] < Cn^r^. 

xG.4„ 

Proof. By theorem 16.71 it follows that there exists constants c,C > such that 

min{cnr2,J^/2} < inf £[x\ < mm{Cn^T\j^} < Cn\^ 

xe^n 

Therefore, if r < J-'/(2c)n^^/^ it follows that min{cnr^, J-'/2} = cnr'^ and the result 
follows. D 
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n 



2_2 



cnr^ < inf f^[x] < Cn^r 




— < inf £'^[x] < jP" 

2 xeAn 



Figure 5. A schematic of regions in the n — T plane in which minimizing deformations 
in An "crossover" to different scahng regimes. In the upper left region minimizing 
deformations in An have larger energy than minimizers over A. the existence of this 
upper bound means that there can be no refinement with decreasing thickness of 
the number of waves for minimizers in A. In the lower left region the minimizing 
deformations scale in energy like that of an isometric immersion. In the region on the 
right the minimizing deformations are either a flat deformation or are close to being 
a flat deformation. It is important to note that this figure is just a schematic drawn 
for a specific values of the constant c in theorem 16.71 Different values of c will change 
the positions of the boundaries between the regions but will not change the qualitative 
form of the figure. 



The previous corollaries 16.81 and 16.91 extend the results of theorem 16.71 and quantify 
different "crossover regimes" in n and r. Specifically, corollary 16.81 gives a critical wave 
number n* > 2 such that energetically there can be no refinement with decreasing 
thickness of the number of waves greater than n*. Corollary 16.91 gives a crossover 
condition for when minimizers transition from being close to a flat deformation to one 
whose elastic energy scales with r like an isometric immersion. In figure [5] we plot a 
schematic of these crossover regimes. 



7. Boundary layer analysis 

The n-periodic isometric immersions constructed in section H] satisfy the FvK equations 
( 1T3ll and (fTSll but not the boundary conditions ( 1161) and (fTTl) . But, in figure H] we see 
that in the vanishing thickness limit numerical minimizers of Sr over An converge to an 
isometric immersion. Therefore, we expect the sheet to be a perturbation of an isometric 
immersion that introduces boundary layers near the edge of the disk and along lines of 
inflection. By theorem 16. 71 the width of these regions will scale with r so that the energy 
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in these regions scales like r^ or smaller. 

Again, let n G {2, 3, . . .}. Define Sn to be the sector in M^ bounded between 6 = 
and 6 = n/n and define Bn = B (1 Sn- Let x G An with corresponding out-of-plane 
displacement rj and potential $ and suppose r] and <l> are perturbations of an n-periodic 
isometric immersion. That is, on Bn assume rj and $ are of the form 



r] = y i^x - cot i^- ) y ) + Tj, 
$ = 1), 



(26) 

(27) 



for some perturbations $ and fj. Furthermore, since x G An we enforce the Dirichlet 
boundary conditions r] = along the lines 6 = and 6 = -. 
The elastic energy of this perturbation is 

2 



^^fxl = 2n (a^) dxdy + 2nT^ / ( Afi - 2 cot (- 

Jb^^ + v ^ ^ Jb„V^-v ^ ^n 

dxdy. 



,"ir\ d'^ri d'^ri 

+4 cot ( - ) ^ + 4^^^ -2[fj,f]] + 2 



nJ dx"^ 



dxdy 



(28) 



If fj and $ extremize the Sr then by the FvK equations (11311171) . and the fact that 6t] = 
along the lines ^ = and 6 = tt/u, it follows that fj and $ satisfy the following boundary 
value problem: 



2(1 + z/) 



4(1 - u] 



.A2$ - 2 cot 



7r\ d 1] d rj 



A^fj + cot - 



n/ dx"^ dxdy 
7r\ d"^^ (92$ 



1 (92$ 5$ 
r d6'^ dr 



nJ dx^ dxdy 



0, 



r=ro,l 



drd6 r d6 



+ [v^v] =0, 

[$,77] = 0, 
= 0, 



1-u 



Arj — n ■ D rj ■ n 



r=ro,l 



' ^^UlM-D^^.t) 



1 — z/ (9n dt 



— CSC 



r=T-o,l 



TT 



n 



r=ro,l 

(z/-2)cos(^)+z/cos(^-2g) 
1-1/ 



2 /7r\ /TT 

- CSC — COS 26 

r \nJ \n 



I dy'^ ' 9x9y / 



^le=0,7r/^ 



9 r] 



9^2 



=0,^ 



6l=0,7r/n 



0, n- 



2r2 cot 



r=ro,l 
521, ^2|,^ 



c)x(9?/ 9x2 



=0,^ 



TT 



n 



(29) 
(30) 

(31) 
(32) 
(33) 

(34) 
(35) 



Remark 7.1. 5*0 far no approximations have taken place. The boundary value problem 
/ f^Plf^) is the Euler-Lagrange equations corresponding to the ansatz l[M\) and [2j 
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Figure 6. The domain for tlie governing equations of the perturbation from an 
isometric immersion is i?„ = {{r,9) : vq < r < 1, < < n/n}. 



Remark 7.2. The boundary value problem (2M3o\) forn = 2 is different from the other 
cases since cot(7r/2) = 0. We will see that for the case n = 2 that there is no need 
to introduce a boundary layer near the lines of inflection. Furthermore, the boundary 
layers near the edges of the annulus will have a different geometry as well. 

We look for approximate solutions of ( 129|) and (!30|) that are linear combinations 
of boundary layer solutions near the interior radius (int), outer radius (out), bottom of 
the sector (bt), and the top of the sector (tp) (see figure [6]). That is we assume 

V = flint + flout + fibt + ^tp and l> = ^int + ^out + ^bt + $tp, (36) 

where each term is found by an appropriate asymptotic expansion. The full configuration 
with domain B can then be obtained by taking odd extensions of t] and even extensions 
of $. 

Remark 7.3. To construct a complete asymptotic solution of the boundary value 
problem ^SMEB^j we would also need to analyse regions where boundary layers overlap. 
In this paper we are interested only in the scaling of the width of the boundary layer and 
do not analyse these overlap regions in this work. 

7.1. Boundary layer near outer radius 

Define the rescaled radius font and functions r^^^^ and ^'^^t by 

font = r'^il - r), flout = r^v'out, ^out = r^Kut 
where — a,/3,7 G M^. To lowest order the stretching and bending energies near the 
outer radius are 



5[x] 
2n 



^27+4« 




?:2 
out 



2n In I- V \ art,, \n 
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+4r^+2° cos(^) (cot (-) + sin(^) 



'^'^out 



+ 2 



dxdy, 



and the compatibility condition is 



1 



Aif^i 



a^^ 



2il + u) 



-T 



d'Vout 



n+^a]^oni _ 2^/3+2a ^^^^g^ {cot{n/n) + sin(^)) 



0. 



out ^ out 

To ensure that the elastic energy is O(r^) and the compatibility condition is non- 
trivial we must have that 

/3 + 2a = 0, 27 + 4a = 2, 7 + 4a = /3 + 2a. 
The solution to these equations gives us the following scaling 

a = -i, /3 = 1, 7 = 2, 
which motivates the asymptotic expansion 



(0) , I (1) , 2 (2) , 



$ 



out 



r'^t + rl^, 



r^$S + .... 



Therefore, to lowest order the equations for the perturbation ( l29l) and ( l30l) become 



54$ 



(^^out 



g^A 



out 



- 4(1 + u) CSC (-) cos{9) cos (-- 9 
+ 4(1 - z/) CSC (-] cos{9) cos (-- 9 



d'^1 



Qf2 



out 



2^(0) 
out 



^2$ 



'^^out 



0. 



Furthermore, the boundary condtions ( 1311) . (1321) and ( l33l) to lowest order become 



9$ 



(0) 



UTout 



0, 



^3j^ 



(0) 
out 



f=0 



drl 



out 



0, 



f=0 



'^'^out 



fout=0 



-2u CSC C-) cos(6') cos (-- 9^ + 2 cot j' 



(37) 
(38) 

(39) 
(40) 



To solve the boundary value problem ( I371H01) we make the following ansatz 
r/(°l = Xl{9)A{r^A-\9)) - B{fout,9), 

^'^^, = Xl{9)C{foutX-\9))-D{f^ 
where A„ : (O, ^) -^ 

/7r\ /TT 

^n{9) = CSC I — ) COs(6') COS 



'out, ' 



is defined by 

The solution to the resulting differential equations containing exponentially decaying 
terms and is given by 



A(f,^,K\9)) 



V 



TT 



V2foutil-i^')^ . l V2fo^t{l-u^)-. 
v/2(l-z.2)'''Pi U9) j'^^i \n{9) 4 
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\n{0)^/2{l-ly^) V / 



^foutil - iy')-A (V2f{l-u^)\ 



cos 



2(1 -z/) "\^ A^^) ) \ \n{e) 

Difouue) = -^^^l^e.p (-V2f{l - .^)iAl(^)) 
A„(6')V2(1 -z/j ^ ^ 

X cos (y2fo«t(l - z/')^a|(^) - ^) . 

A^(^)A(fo„iA~^(^)) and \1^{9)C {vouAn^ (0)) can be interpreted as terms that alone 
would reduce the magnitude of the Gaussian curvature in a thin boundary layer while 
B{fout^d) and D{rout,(^) would alone locally reduce the mean curvature in a separate 
boundary layer with a different geometry. In the overlap of these two boundary layers the 
bending energy is reduced through the combination of these two effects. By expressing 
?7qJj and $J,J^ in terms of the actual radius p = Rr = Rr^/'^f we have the following 
results for the scaling of the width of these boundary layers: 

(i) The width of the boundary layer in which the Gaussian curvature is significantly 
reduced satisfies the following scaling 

Wv\i\i{9) p=R ~ t2 |/^q|~4 csc(7r/'ri) cos(^) cos(7r/n — 9). 

(ii) For n > 3, the width of the boundary layer in which the mean curvature is 
significantly reduced satisfies the following scaling 

width(6')p=^ ~ tali^oT* V'sin(7r/n) sec(6') sec(7r/r 



'n 



7.2. Boundary layer near interior radius 

The boundary layer near r = ro is completely analogous to the one near r = 1. Define 
the rescaled radius fint by 

Tint = r"2(r-ro) 
and consider the asymptotic expansion 

(0) , i (1) , 2 (2) , 
Vint = m^t + '^'mA + ^ Vint + • • • 

^int — T <P.„^ + r2 <Pj„^ + T <Pj„4 + 

This yields identical governing and boundary equations as the outer radius and thus 
77J2 = Xl{9)A{n^t, X-\9)) - B{hnt, 9) 

^^t = ^liO)A{r,r.t), X-\9) - D(f,„„ 9), 
where A,B,C,D are defined as in the previous subsection. 
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7.3. Boundary layer near the bottom of the sector 



28 



From the observations in the remark 17.21 we will assume in this section that n G 
{3,4, . . .}. Define the rescaled coordinate ybt and functions t^^^ and $'f,^ by 

ybt = T'^y, Vbt = r^v'bt, ^bt = T^^'i^t 

where — a,/3,7 G M"*". To lowest order the stretching and bending energies near y = 
are 

2 



5[x] 
2n 



T 



27+4a 



1 



3^r 



92$ 



bt 



T 



m^] 



2n 



T 



B ^ + 1^ 
1 , 



dy, 



-.1 

bt 



dxdy, 



Bn 



l-v 



./3+2a 



dylt 



and the compatibility condition is 



-r 



7+4o 



a^$ 



bt o_/3+a ^ Vbt 



-2t^ 



2 cot ( - 



0. 



dxdy. 



2(1 + z/)- dyl, -■ dxdyu 

Therefore, the scaling that ensures the elastic energy is O(r^) and the compatibility 
equation is non-trivial is 



a 



/3 



2 ^-5 



3) K 3' ' 3' 

which is a different scaling then the one near the edges of the annulus. This scaling 
motivates the asymptotic expansion 



Vbt -- 

^bt 



(1) , ^ (2) , 



2 (0) 
T'Vbt 



Consequently, to lowest order (|29|) and (!30|) become 



dybt 

d 
dm 

Furthermore, by 



dylt 



bt _^l + ^)Ht 



dx 






dy-i 



r,3 

bt 



dx 



0, 



0. 



(41) 
(42) 



and 035p the boundary conditions are 



2^(0) 



52$ 



bt 



dx"^ 



2^(0) 



52$ 



bt 



ybt=o 



(0) 



dxdyu 



ybt=o 



2/64=0 



5^i 



r,2 



?;6t=0 



/vr 

2 cot - 

Vn 



To solve this boundary value problem we make the following ansatz 

Vh = Mx) + ylt cot {l) fix, ybt), 
^'^^=ilj,{x)+ylcot{^)g{x,ybt), 



(43) 
(44) 

(45) 
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which transforms equations ( HTj) and (H2]) into 

i: fe& + 6y..& + 6^ - 4(1 + u)yl'^) = 0, 



dytt 



'■dvl 



■ m 



dytt 



i- y!M + 6y,,^ + 6^ + m-u)yl 



btdx 



0. 



29 



(46) 



dytt V ^^% ^Vbt ^ybt 

If we integrate with respect to ybt and make the similarity transformation z = y^^/x, we 
have the following ordinary differential equations 

27.^& + 108.& + 60g + 4(l + .).^ = 0, 
27.20 + 108.S + 60£- 4(1 -.).g = 0. 

Solving for -^ we obtain the following single differential equation 



,, df /90d/ 
4(1 -z.2).-^+ --^ + 2430 



dV .... dV 3645 2dV 729 .d^/ 

1-4 + 4455.-4 + -2^4 + ^^ ^ 

d^2 d^3 2 d.4 4 d.5 



, , , „ , ~ ,0. (47) 

dz'Vzdz j~2' j^.s' o j„4' /fj^rs/ vv 

The general solution to equation ( 1471) that does not contain exponentially growing 



terms is 



/(^) 



ci."tKeri (a ■ 3~^ ■ (1 - u^)^^/z) + caz'^Keii f4 ■ 3~^ • (1 - u'^)^^/z 



dz 



+ C3, 



where Keri/3 and Keii/3 denote Kelvin functions of the second kind and ci, C2 and C3 
are arbitrary constants. Therefore, we have that 



(0) 

Vbt 



ciz 2Keri (4-3 2 ■ (1 



^i{x)+yltCoi{'n/n) 

+C2z"^Keii (a ■ 3"5 . (1 _ z/^)!^ 



y )^\^z 



dz + C3 



$ 



if =^2(x)+yMC0t(7r/n) 



(\Iy^[ h^""^*"'^ (^ ■ ^"' ■ (^ " '"'^'^ 



-ci. aKeii [4-32 



u l^x/z 



dz + C4 



where C4 is an arbitrary constant. 
Now, near y = 0, we have that 



2--32r I - (1 - Z/2)-T5 (C2 - Ci) — 

'^/ Vbt 



Vu =V'i(a;)+ywCot (^^j 

-2-^3^r ( -M (1 - z/')i^ (ci (1 + Vs) + C2 (1 - Vs)) ?^ + C3 + . . . 



$ 



(f = V]>2(x)+y2^COt(^ 



17 3„ / 1 

-2'-32r( - 



u^)-^^ 



2 

X3 



(1 - U) (Ci + C2) -32" 



13 1 „ 
-2^-33! 



+ C4 



I) (^ 



z/ 



0^(1 



v-1 



-l + y3)ci+ (1 + 73)02) — 



1 

X3 



Vbt 
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Furthermore, 

lim / ^-iKeri (a ■ 3^5(1 - u'^)^^] dz = -J2 + Vs, 
2—^00 J 3 v / 2 

lim / ^"tKeii [4 ■ 3"i(l - v'^f^^^ dz = -\J2 - Vs. 
Therefore, to satisfy the boundary conditions (H3!l and (1^ we must have that 



3_/-, X /_N _,„,!. 






•^1 ~ — ^: — rn~5 ^2 — -3- — — T— , l^c 



V6(l-i^2)i7r' 3^(l-i/2)^7r' 



C3 = 1, C4 



1+;^ 



3(1-1/)- 

By expressing 77^^ and $^j in terms of the actual Cartesian coordinates u = Rx 
and V = Ry and approximating p and 9 near y = by p ~ m and 6* ~ v/u, it follows 
that the width of the boundary layer in this region scales like 



" 6 



width(p)^=o ~ t'^ip^lKol 

This boundary layer can be interpreted as a region in which the mean curvature is 
locally reduced while the change in energy contributed from the Gaussian curvature of 
the perturbation is of the order 0(r^/^). 

7.4- Boundary layer near the top of the sector 

To construct the boundary layer near 6 = n/n we can simply rotate the asymptotic 
expansion for fj^t and "I'm through the angle 6 = Tc/n and then evenly reflect about the 
line 6 = ix/n. That is, in polar coordinates we set 

Vtpir, 0) = fju (r, Tx/n - 6) and l>tp(r, 6) = ^^ (r, n/n - 9) . 

8. Discussion 

In this paper we studied the convergence and scaling with r for minimizers of the Foppl- 
von Karman energy for non-Euclidean plates with constant negative Gaussian curvature 
Kq. Specifically, to obtain a better understanding of the experimental results in [20] , 
we focused on annular domains and deformations with a periodic profile of n waves. 

The first main result of this work is theorem 16. 7[ which gives rigorous upper and 
lower bounds for the elastic energy of minimizing deformations. This theorem rigorously 
proves that for all thickness values the 2-wave saddle shape is energetically preferred 
over deformations with a higher number of waves. Specifically, this theorem proves that 
in the FvK approximation there can be no refinement of the number of waves with 
decreasing thickness in contrast to what is observed experimentally. 

The second main result is the scaling with r of the width of boundary layers in 
which stretching energy is concentrated. In these localized regions of stretching, the 
total elastic energy of an isometric immersion is lowered by allowing some stretching 
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to reduce the bending energy. Near the edge of the annulus the contributions to the 
bending energy coming from the Gaussian and mean curvatures is reduced in overlapping 
boundary layers that scale with the thickness like t^^'^\Ko\~^^^. Furthermore, for n > 3 
boundary layers form along the lines of inflection in which the mean curvature of an 
isometric immersion is reduced by removing a jump discontinuity in the curvature along 
the azimuthal direction. The width of these boundary layers scales with the radius and 
thickness like p^/^t^/^\Ko\-^/^. 

It is important to note that the FvK elastic energy is valid in an asymptotic limit 
in which the metric becomes increasingly flat with decreasing thickness, that is e ~ t/R. 
In this paper we took e fixed and t decreasing and in particular we showed that the 
energy of minimizers scales like r^. Consequently, when r ^ 1 it is perhaps more 
appropriate to consider the Kirchhoff model. But, the Kirchhoff model is rigid in the 
sense that there is no competition between stretching and bending energies and thus 
there is no possibility of the Kirchhoff model alone explaining the refinement of the 
number of waves with decreasing thickness. 

It should also be pointed out that by the results of Hilbert [31], Holmgren [35] 
and Amsler [36] concerning that non-existence of analytic isometric immersions of the 
hyperbolic plane into M^, that -Eki diverges in R for fixed local isometric immersions. 
This is a geometric feature of non- Euclidean plates that is not captured in the linearized 
geometry of the FvK model. Indeed it was conjectured in [19] that 

max 1 1 A^i I , \k2\} > exp ( — - ) , 

where ki,k2 are the principal curvatures of the surface x(P). This is in contrast with 
the FvK model in which \D'^r]\ does not grow with the size of the domain. 

Furthermore, it was shown in [19] that the n— periodic isometric immersions in 
the FvK ansatz approximate exact isometric immersions with a periodic profile. The 
principal curvatures - and thus the bending energy ii^Ki ^ of these exact isometric 
immersions diverges at a finite radius i?„ that scales with n like -R„ ~ log(n) [19]. This 
gives a geometric mechanism for the refinement of the number of waves with increasing 
radius of the disk but it does not explain the experimentally observed refinement with 
decreasing thickness. 

From this work and [19] it is clear that the FvK and Kirchhoff models of non- 
Euclidean elasticity cannot completely explain the periodic shapes in [20] . The periodic 
shapes we have constructed in this paper are qualitatively similar to the experimental 
shapes but are only local minimizers of the FvK energy. A major goal of future 
works is to connect the existence of these local minimizers to the observed patterns 
in experiments. Below we outline two avenues of future research that could shed light 
on these issues. 

First, the exact isometric immersions with periodic profile constructed in [19] are 
similar to the n-periodic isometric immersions in the FvK but have highly localized 
regions of bending energy near the edge of the disk and along the lines of inflection. 
Therefore, in these regions, as in the boundary layers in the FvK ansatz, it is conceivable 
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that it would be energetically favourable to allow some stretching to reduce this localized 
bending energy. These observations illustrate the multiple scale behaviour of this 
problem - namely the scales t/R and ^/\Ko\R - and it may be more appropriate to 
consider a combination of different reduced theories in various regions of the domain. A 
hierarchy of such reduced theories has recently been conjectured by Lewicka, Pakzad, 
and Mahadevan [15] and further research in this direction may explain the complex 
morphologies of non-Euclidean plates. 

Second, the periodic shapes in swelling hydrogels are the result of dynamical 
processes. It may be more appropriate to model this type of differential growth 
dynamically, perhaps as a gradient flow of the elastic energy. The pattern could then be 
selected for dynamical reasons and not by global minimization of an energy functional. 
This might explain why local but not global extrema for the energy functional seem to 
describe the observed patterns. 
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